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Abstract 

We discuss the phenomenological viability of string multi-Higgs doublet models, namely a 
scenario of heterotic Z3 orbifolds with two Wilson lines, which naturally predicts three supersym¬ 
metric families of matter and Higgs fields. We study the orbifold parameter space, and discuss 
the compatibility of the predicted Yukawa couplings with current experimental data. We address 
the implications of tree-level flavour changing neutral processes in constraining the Higgs sector 
of the model, finding that viable scenarios can be obtained for a reasonably light Higgs spectrum. 
We also take into account the tree-level contributions to indirect CP violation, showing that the 
experimental value of ek can be accommodated in the present framework. 


1 Introduction 

The understanding of the observed pattern of quark and lepton masses and mixings remains as one 
of the most important open questions in particle physics. From experiment, we believe that Nature 
contains three families of quarks and leptons, with peculiar mass hierarchies. Moreover, there is firm 
evidence that the flavour structure in both quark and lepton sectors is far from trivial, as exhibited 
by the current bounds on the quark [1] and lepton [2-4] mixing matrices. 

The standard model of strong and electroweak interactions (SM) fails in explaining some important 
issues such as the fermion flavour structure or the number of fermion families we encounter in Nature. 
Moreover, the mechanism of mass generation for quarks and leptons is still unconfirmed, since the 
Higgs boson is yet to be discovered in a collider. 

Other high-energy motivated theories, such as supersymmetry (SUSY), supergravity (SUGRA), or 
grand unified theories (GUTs) may repair some shortcomings of the SM (as for instance the hierarchy 
problem, the existence of particles with distinct spin, or the unification of gauge interactions), but 
they still fail in providing a clear understanding of the nature of masses, mixings and number of 
families. In this sense, a crucial ingredient to relate theory and observation is the precise knowledge 
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of how fermions and Higgs scalars interact, in other words, the Yukawa couplings of the fundamental 
theory. 

String theory is the only candidate to unify all known interactions (strong, electroweak and grav¬ 
itational) in a consistent way, and therefore it must necessarily contain the SM as its low-energy 
limit. In this sense, string theory must provide an answer to the above mentioned questions. A 
very interesting method to obtain a four dimensional effective theory is the compactification of the 
Eg x Eg heterotic string [5] on six-dimensional orbifolds [6], and this has proved to be a very successful 
attempt at finding the superstring standard model [7-28] (other interesting attempts at model build¬ 
ing using Calabi-Yau spaces [29], fermionic constructions [30,31], and heterotic M-theory [32,33], 
can be found in Refs. [34-36], [37-42], and [43], respectively). As it was shown in [8,10], the use 
of two Wilson lines [6, 7] on the torus defining a symmetric Zg orbifolcl can give rise to SUSY mod¬ 
els with SU{2>) x SU( 2) x U(l) n gauge group and three families of chiral particles with the correct 
SU (3) x SU (2) quantum numbers. These models present very attractive features from a phenomeno¬ 
logical point of view. One of the U( l)s of the extended gauge group is in general anomalous, and it 
can induce a Fayet-Iliopoulos (FI) ZUterm [44-47] that would break SUSY at very high energies (FI 
scale ~ 0(1O 16 - 17 GeV)). To preserve SUSY, some fields will develop a vacuum expectation value 
(VEV) to cancel the undesirable D-term. The FI mechanism allows to break the gauge group down 
to SU( 3) c x SU(2)l x U{1)y and obtain the mass spectrum of the minimal supersymmetric standard 
model (MSSM), plus some exotic matter, as extra singlets, doublets or vector-like triplets, depending 
on the model, as shown in Refs. [15,16] and [11]. 

Orbifold compactifications have other remarkable properties. For instance, they provide a geome¬ 
tric mechanism to generate the mass hierarchy for quarks and leptons [48-52] through renormalisable 
Yukawa couplings. Z n orbifolds have twisted fields which are attached to the orbifold fixed points. 
Fields at different fixed points may communicate with each other only by world sheet instantons. The 
resulting renormalisable Yukawa couplings can be explicitly computed [48,49,53-56] and they receive 
exponential suppression factors that depend on the distance between the fixed points to which the 
relevant fields are attached. These distances can be varied by giving different VEVs to the T-moduli 
associated with the size and the shape of the orbifold. 

However, the major problem that one encounters when trying to obtain models with entirely 
renormalisable Yukawas lies at the phenomenological level, and is deeply related to obtaining the 
correct quark mixing. Summarising the analyses of Refs. [51,52], for prime orbifolds the space group 
selection rules and the need for a fermion hierarchy forces the fermion mass matrices to be diagonal at 
the renormalisable level. Thus, in these cases, the Cabibbo-Kobayashi-Maskawa (CKM) parameters 
must arise at the non-renormalisable level. For analyses of non-prime orbifolds see Refs. [51,52,57,58]. 

For example, since the FI breaking generates VEVs for fields of order (xj) ~ 10 16 ~ 17 GeV, if 
one has terms in the superpotential of the type \i • ■ ■ Xm £ £ £, these would produce couplings of 
order ( Mpi/Mp) m . Therefore, depending on m, different values for the couplings might be generated. 
Obviously, the presence of these couplings is very model-dependent and introduces a high degree of 
uncertainty in the computation. However, it is important to remark that having the latter couplings 
is not always allowed in string constructions. First of all, they must be gauge-invariant, something 
that is not easy to achieve, due to the large number of U( 1) charges which are associated to the 
particles in these models. Even if the couplings fulfil this condition, this does not mean that they are 
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automatically allowed. They must still fulfil the so-called “stringy” selection rules. For example in the 
SU( 3) x SU( 2) x U(1)y x SO(10)hidden model of Ref. [15], where renormalisable couplings are present, 
only a small number of non-renormalisable terms are allowed by gauge invariance. Nevertheless, even 
the latter terms turn out to be forbidden by string selection rules. 

Clearly, purely renormalisable Yukawa couplings are preferable, because, in general, due to the 
arbitrariness of the VEVs of the fields entering the non-renormalisable couplings, the predictivity is 
lost. Furthermore, as discussed above, higher-order operators such as those induced by the FI breaking 
are very model-dependent. One possibility of avoiding the necessity of these non-renormalisable 
couplings is to relax the requirement of a minimal matter content (with just two Higgs doublets) in a 
Z 3 orbifold with two Wilson lines. Since these models naturally contain three families of everything, 
including Higgses, additional Yukawa couplings will be present, with the possibility of leading to 
realistic fermion masses and mixings, entirely at the renormalisable level (with a key role being played 
by the FI breaking) [25]. In addition, and given the existence of three families of quarks and leptons, 
having also three families of Higgses renders these models very aesthetic. In fact, let us recall that 
experimental data imposes no constraints on the number of Higgs families. Moreover, this non-minimal 
Higgs content, provided that the extra doublets are light enough to be present at low-energies, also 
favours the unification of gauge couplings in heterotic string constructions. Due to the FI scale, the 
gauge couplings may unify at the string scale (« (/gut x 5.27 10 17 GeV) [24]. 

Thus, this class of string compactifications is one of the scenarios where one can obtain a SM/MSSM 
compatible low-energy theory, albeit with an extended Higgs sector. Furthermore it offers a solution 
to the flavour problem of the SM and MSSM, since the structure of the Yukawa couplings is completely 
derived from the geometry of the high-energy string construction. Given the increasing experimental 
accuracy, accommodating the data on quark masses and mixings is not straightforward. In this work, 
we propose to investigate in detail whether or not it is possible to obtain Z 3 orbifold configurations 
that successfully reproduce the observed flavour pattern in Nature. In this sense, having additional 
Yukawa couplings presents several advantages, as for example a greater flexibility when fitting the 
data from the quark masses and mixings. 

On the other hand, when working in a multi-Higgs context, we should also take into account the 
potential appearance of flavour-changing neutral currents (FCNCs) at the tree level, which could con¬ 
tribute to a wide variety of Higgs decays and interactions with other particles [59-76]. Generally, the 
most stringent limit to these flavour-changing processes is assumed to come from the mass difference 
of the long- and short-lived neutral kaons, A ttik = — m K s ■ A possible way to overcome this 

problem is by imposing that the Higgs spectrum is heavy enough to suppress the undesired contribu¬ 
tions to the neutral meson mass differences. As we will see, the Yukawa couplings of this Z 3 scenario 
exhibit a strongly hierarchical structure, and this property is instrumental in circumventing the FCNC 
problem without the need for an excessively heavy Higgs sector. 

Z 3 orbifolds are also very attractive when addressing the lepton sector, and in fact offer an ap¬ 
pealing scenario to study the problem of neutrino masses (predicting naturally small Dirac masses, 
and thus a low see-saw scale). We postpone this analysis to a forthcoming work [77]. 

This work is organised as follows. In Section 2, we describe the main properties of the Yukawa 
couplings in Z 3 orbifold models. We study the relations between the several orbifold parameters 
induced from the quark mass hierarchy and from electroweak symmetry breaking. Section 3 is devoted 
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Figure 1: Two dimensional sublattices (i = 1, 3, 5) of the Zj, orbifold, symbolically denoting the fixed 
points as (o, •, x). 


to a brief overview of the extended Higgs sector. In Section 4 we present the contributions of neutral 
Higgs exchange to tree-level FCNCs. The numerical analyses of the orbifold parameter space and 
FCNCs in association with specific Higgs textures is given in Section 5, where we also address the 
possibility of new contributions to indirect CP violation. Finally, we summarise our results in Section 6. 


2 Yukawa couplings in Z% orbifold models 

In this Section we review some of the most relevant features of the geometrical construction of the 
Z^ orbifold leading to the computation of the quark mass matrices. We study the correlations of the 
orbifold parameters arising from the quark mass hierarchy and electroweak (EW) symmetry breaking 
conditions, and derive useful relations which play an important role in constraining the parameter 
space. 

2.1 Z$ orbifold: a brief review 

The Z-i orbifold is constructed by dividing R 6 by the [S'!/(3)] 3 root lattice modded by the point group 
(P) with generator 9 , where the action of 9 on the lattice basis is 9ei = ej+i, 9ei+\ = —(e, + e,+i), 
with i = 1,3,5. The two-dimensional sublattices associated to [Sf7(3)] 3 are presented in Fig. 1. 
Invariance under the point group reduces the orbifold deformation parameters to nine: the three radii 
of the sublattices and the six angles between complex planes. The latter parameters correspond to 
the VEVs of nine singlet fields appearing in the spectrum of the untwisted sector, and which have 
perturbatively flat potentials. These so-called moduli fields are usually denoted by T. 

In orbifold constructions, twisted strings appear attached to fixed points under the point group. 
In the case of the Z 3 orbifold there are 27 fixed points under P, and therefore 27 twisted sectors. We 
will denote the three fixed points of each two-dimensional lattice as in Fig. 1. In the Z 3 orbifold, the 
general form of the Yukawa couplings between the twisted fields is given by a Jacobi theta function, 
and their expressions can be found, for example, in the Appendix of Ref. [56]. The Yukawa couplings 
contain suppression factors that depend on the relative positions of the fixed points to which the 
the fields involved in the coupling are attached, and on the size and shape of the orbifold (i.e. the 
deformation parameters). Let us first study the situation before taking into account the effect of 
the FI breaking. Let us suppose that the two non-vanishing Wilson lines correspond to the first and 
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second sublattices. Then, the 27 twisted sectors come in nine sets with three equivalent sectors in each 
one. The three generations of matter (including Higgses) correspond to changing the third sublattice 
component (o,*, x) of the fixed point, while keeping the other two fixed. Consider for example the 
following assignments of observable matter to fixed point components in the first two sublattices, 


Q <-> o o u c <-> o o 

L • e c <-> • x 

H u <-> o o H d 


d c <-> x o 
v c <-> x x 


In this case, the up- and clown-quark mass matrices are given by [25]: 


where 


A u 


M u = gNA u , M d = g N si A d , 


/ 

W 2 

we £5 

W 4£5 ^ 

, A d = 

( 

W\ 

we £5 

w 3 £5 


we £5 

W4 

W 2 £5 


we £5 

W3 

Wl £5 

V 

U> 4 £ 5 

W2 £5 

We ) 


V 

1^3 £5 

Wl £5 

w 5 


( 1 ) 

( 2 ) 

(3) 


g is the gauge coupling constant, and N is related to the volume of the Z 3 lattice unit cell such that 
gN ss 1. In the above matrices Wi denote the VEVs of the neutral components of the six Higgs 
doublet fields 1 . Since we are assuming an orthogonal lattice, i.e. with the six angles equal to zero, 
only the diagonal moduli (T)), which are related with the radii of the three sublattices 2 , contribute to 
the Yukawa couplings, through the suppression factors £* 


£i ~ 3e a £ Ti i,j = 1,3,5. 


(4) 


2.2 Quark mass matrices and Yukawa couplings after the Fayet-Iliopoulos 
breaking 

As mentioned in the Introduction, the anomalous U( 1) of the extended SU( 3) x SU{2) x U(l) n gauge 
group generates a Fayet-Iliopoulos H-ternr which could in principle break SUSY at energies close 
to the string scale. This term can be cancelled when scalar fields (C,), which are singlets under 
SU(3) x SU(2), develop large VEVs (10 16-17 GeV). The VEVs of these fields (c*), have several 
important effects. Firstly, they break the original SU(3) x SU{2) x U( 1)" gauge group down to 
the (MS)SM SU( 3) x SU( 2) x U( 1). Secondly, they induce very large effective mass terms for many 
particles (vector-like triplets and doublets, as well as singlets), which thus decouple from the low-energy 
theory. Even so, the SM matter remains massless, surviving as the zero mass mode of combinations 
with the other (massive) states. All these effects modify the mass matrices of the low-energy effective 

1 For convenience, we adopt a distinct notation from the one originally used in Ref. [25]. Instead of denoting the 
Higgs fields (and the corresponding VEVs) Hf’ u (vf ,u ), i = 1,2,3, we use Hi (), i = 1,... ,6, with odd (even) cases 
referring to down- (up-) quark coupling Higgses. 

2 Actually, in Ref. [25] it is argued that it is possible to fit the entire fermion data by using only one degenerate 
radius. In our approach, and given the increasing accuracy in the experimental data, we allow for the more general case 
of three radii (and thus three moduli). 
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theory (see Eq. (2)), which, for the example studied in [25], are now given by 3 


M u =gNa u A u B u , 

M d = g Nei of A d B d °, 

where A u,d are the quark mass matrices prior to FI breaking (see Eq. (3)), af is given by 

r f 

f C 9 

a J = 


c {| 2 + |c {| 2 


with / = u c ,d c , and B? is the diagonal matrix defined as 

B f = diag(/3 / e 5 , l,a f /e 5 ). 

Finally 


a f = £51 


I 2 + | c f J 2 


\c{s 5 \ 2 + \c f 2 \ 2 


V = 


\ 4\ 2 + \ 4\ 2 


of 12 


+ |C2£5 


f Cr 12 


( 5 ) 


( 6 ) 


(7) 


( 8 ) 


In the above, cj are derived from the VEVs of the heavy fields responsible for the FI breaking as 


c{ = e ' {f) c! 


of = F "U) r f 
to - C Eo 


— c j ^2 — c L 2 5 

where in each case e' and e" can take any of the following values: 

e ', e" = 1 , ei, £ 3 , £1 £3 . 


(9) 


( 10 ) 


Let us also stress that one should not take a/, pf , £5 and a.f as independent parameters. In fact, 
Eqs. ( 6 , 8 ) imply that 


a? = 


(l - a /2 ) 


1/2 


£5 


a f (1 -£§) 1/2 

so that for given values of £5 and af, fif is fixed as 

P f = 


= 1 - 




1/2 


U-£§) 1/2 ’ 




(ii) 


( 12 ) 


Eqs. (5,7) become more transparent when the terms that encode the flavour structure are explicitly 
displayed: 


• (13) 



/ 

w 2 £5 P uC 

W 6 £5 

u>4 a“° ^ 


/ 

wi £5 P dC 

w 5 e 5 

W 3 Oi dC \ 

A u B uC = 


w 6 e\ (3 uC 

W4 

IC 2 a u 

, A d B dC = 


w 5 £5 P dC 

w 3 

w\ a dC | 


V 

w A e\ P uC 

ui 2 £5 

w 6 a uC /e 5 y 


V 

w 3 el P dC 

Wl £5 

w 5 a dC /e 5 ) 


3 Note that although the Ci are in general complex VEVs, they only introduce a global and therefore unphysical phase 
in the mass matrix. More complicated examples would in principle give rise to a contribution to the CP phase [25]. 
This mechanism to generate the CP phase through the VEVs of the fields cancelling the FI D-term was used first, in 
the context of non-renormalisable couplings, in Ref. [51]. For a recent analysis, see Ref. [78]. 
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Given that the mass matrices are related to the Yukawa couplings as 


M u = WiY i U ’ M d= H w i y i i 

i=2,4,6 2=1,3,5 


( 14 ) 


the structure of the Yukawa couplings is easily derived from Eq. (13). For the down sector, the latter 
read: 



/ 

£5 (3 dC 

0 

0 ^ 

{ ° 

0 

a dC \ 

Y d = gNe ia dC 


0 

0 

a dC , Y d = gN £l a dC 

0 

1 

0 


V 

0 

£5 

0 ) 

V 4 P dC 

0 

0 ) 


Yr“ 


gNeia dC 


( 0 

4 P dC 

V o 


£5 0 

0 0 

0 a d °/ e 5 


(15) 


The Yukawa couplings for the up-type quarks can be also obtained by doing the appropriate replace¬ 
ments: (eia d °) —> a“ c and a d °, j3 d ° —> a“°,/3“ c . 

Expanding the eigenvalues of the quark mass matrices up to leading order in e 5 , one can derive 
the following relation 4 for the Higgs VEVs in terms of the quark masses 5 [25] 


down-quarks: 
up-quarks: 


{w 1 ,w 3 ,w 5 j (gN e i a d ' 
{w 2 ,w 4: ,w 6 } (gN a u ) 



,m s , 


m b £ 5 \ 
a dc / 


, m c 


m t £5 


a u 


(16) 


The most striking effect of the FI breaking is that it enables the reconciliation of the Yukawa 
couplings predicted by this scenario with experiment. In particular, and as we will see in Section 5.1, 
the quark spectra and a successful CKM matrix can now be accommodated. 


2.3 EW symmetry breaking and the orbifold parameter space 


In addition to the hierarchy constraint imposed by the observed pattern of quark masses, the VEVs 
must further comply with other constraints as those arising from EW symmetry breaking (EWSB): 

w\ + w\ + uig + wl + w\ + Wq = 2 Mf/ (g 2 + g' 2 ) « (174 GeV) 2 . (17) 


In particular, we have that 


1 


1 


(: gNa uC ) 2 

1 

(, gNeia dc ) 2 


(£ 5 /3 “ c ) 2 

1 

d c \ 2 


m, 


r m? 

+ ^1 — 

m c 


m„ 


(£ 5 /?' 


m d +e b 5 - 


■ mi 


mz 


\ a u / 

(m b e 5 y 

\ a dc ) 


(174 GeV) 2 


(18) 


4 Regarding quark mixing, it is also possible to obtain analytical expressions (up to second order in £ 5 ) for the several 
CKM matrix elements, as done in Ref. [25]. 

5 Notice that there is a misprint in these equations in Ref. [25], where in the corresponding version of Eq. (16) the 
factor appeared as £ 5 . 
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We notice that the above condition can always be fulfilled since the quark Yukawa matrix prefactors, 
£1 and gN, have not yet been used. At this point, let us introduce a generalised definition for tan/3: 


Vd + U!§ + w% 

Using Eq. (16), Eq. (19) can be rewritten as 


(19) 


tan/3 = £i — 
a u 




( m «+ e iS) + m c+( 


\ ( es /? 1 ^) 2 ( md + £ i^r) + W s + (l^ L ) 


( 20 ) 


From the above equation it becomes manifest that by considering a given value for tan /3 we are 
implicitly defining eq, for fixed values of £5 and oJ. This in turn implies that according to Eq. (18), 
gN is in fact a function of tan/3, £5 and a/, and its value, gN ns 1, suffers tiny fluctuations (of 
order 1% - 10%) in order to accommodate the correct EWSB. The latter statements become more 
transparent noticing that by bringing together Eqs. (18) and (20), one can derive useful relations that 
allow to express gN and £1 as a function of the quark masses and orbifolcl parameters for a given 
value 6 of tan /3: 


gN 


1 (l + tan 2 /3 ) 1 / 2 
a“ c tan /3 



174 GeV 


£igN = — (! + tan 2 0 ) 


1/2 


\j (£5/3 ^) 2 ( TOd + + m s + (1^) 


174 GeV 


( 21 ) 


The first equality of Eq. (21) provides a clear insight to understanding the smallness of the fluctuations 
of gN. Assuming the limit where ,£5 <C 1, a“° ~ £ 5 / 0 “°, so that gN « m t /(174 GeV). 

It is also important to comment on the relative size of the VEVs 61 and 62 . From the definition of 
at (Eq. ( 6 )) we can derive an additional relation 


I c{ I 


£"(/) 



( 22 ) 


where we have used the definitions of Eqs. (9,10). If, for example, one assumes the VEVs to be of the 
same order of magnitude, i.e. Ci ~ C 2 , then one should further ensure that 


£"(/) 

£'(/) 



~ 1 . 


(23) 


To conclude this Section, let us make a few remarks regarding two topics so far not discussed. 
Firstly, and since it is well known that the CP symmetry is not conserved in nature, it is important to 
comment on the sources of CP violation present in this class of models. The Yukawa couplings have 

6 Whenever referring to a parameter whose value was estimated using the EWSB conditions, and which is a function 
of tan /3, we will use the designation “EWSB fit”. 



been defined through real quantities, so that no physical phase appears via the CKM mechanism. 
However, this need not be the most general scenario. Dismissing for the present time the possibility 
of spontaneous CP violation, associated with non-trivial phases of the Higgs VEVs, there still remains 
another source of CP violation, in addition to the one already mentioned in footnote 3. Should the 
VEV of the moduli field have a phase, then CP (which is a gauge symmetry of the model) would be 
spontaneously broken at very high energies. The phases would be fed into (thus also appearing in 
a/), and would be present in the Yukawa couplings. Therefore, in the low-energy theory, CP would 
be explicitly violated via the usual CKM mechanism [79,80]. 

It is also relevant to mention the effect of the renormalisation group equations (RGE) on the mass 
matrices presented in this Section. The flavour structure of Eqs. (5,13) is associated with a mechanism 
taking place at a very high energy scale. However, and given the clearly hierarchical structure of the 
quark mass matrices, one does not expect that RGE running will significantly affect the predictions 
of the model. 


3 The extended Higgs sector 


As mentioned in the previous sections, in this class of orbifold models, one has replication of families 
in the Higgs sector. By construction, this scenario contains three generations of SU (2) Higgs doublet 
superfields, with hypercharge —1/2 and +1/2, respectively coupling to down- and up- type quarks. 


H 


1(3.5) 


h° 

([1(3,5) 

^1(3,5) 


H■ 


2(4,6) 


h + 

"■2(4,6) 

to 

z(4,6) 


(24) 


In Ref. [76], we have studied the general case of SUSY models with Higgs family replication, as is the 
present case. Hence, we will just summarise here some important features which are relevant for the 
present analysis. We assume the most general form of the superpotential, which is given by 


W = Q (Yf H, + Y d H 3 + Y d H b )D c + L (Yj e Hi + Y 3 e H 3 + Y 5 e H 5 )E c 
+ Q {Y?H 2 + Y?Hi + Y“H 6 )U c + Mi 2 ^ii ?2 + M14M4 + rneMe 
+ H32H 3 H 2 + n 3 4H 3 Hi + /j, 3 qH 3 Hq + h$ 2 H 3 H 2 + , ( 25 ) 


where Q and L denote the quark and lepton SU(2)l doublet superfields, U c and D c are quark singlets, 
and E c the lepton singlet. The Yukawa matrices associated with each Higgs superfield, l/ 9 , have been 
already defined in Eq. (15). In what follows, we take the /i,j as effective parameters. (Notice that in 
this context the Giudice-Masiero mechanism [81] to generate the fi -term through the Kahler potential 
is not available for prime orbifolds such as the Z 3 orbifold [82,83].) 

The scalar potential receives the usual contributions from D -, F- and SUSY soft-breaking terms, 
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which we write below, using for simplicity doublet components. 


vf = E E & 


i,j= 1,3,5 
1= 2,4,6 


i=l,3,5 
fc,Z=2,4,6 




a=l 


' 6 

2 

n' 2 

" 6 

Y h 'I t " h ' 

+ T 

E(-^ i^i 2 


Koft = E + Y, (mDuHtm- Y [{Bi*)i j H i H j +H.c.]. (26) 

2,j=1,3,5 k,l— 2,4,6 


i=l,3,5 
J=2,4,6 


After electroweak symmetry breaking, the neutral components of the six Higgs doublets develop 
VEVs, which we assume to be real, 


(^1(3,5)) “ w l(3,5) ! (*-2(4,6)) ~ u, 2(4,6) > 


(27) 


and as usual one can write 


+ 71 (<T * 


- i<Pi) ■ 


(28) 


For the purpose of minimising the Higgs potential and computing the tree-level Higgs mass matrices, 
it proves more convenient to work in the so-called “Higgs basis” [59,65], where only two of the rotated 
Helds develop VEVs: 


4v. Pij hj , 

(</>?) = v d , (0°) = V u ■ (29) 

By construction (cf. Eq. (17)), the new VEVs must satisfy 

v 2 u + v 2 d « (174 GeV) 2 , (30) 

and we can now define tan/3 (see Eq.(19)) in the standard way, 

tan/3 = — . (31) 

v d 

In the new basis, the free parameters at the EW scale are which has dimensions mass 2 , and 

tan /3 (for a detailed discussion of the Higgs basis, including the definition of the new parameters and 
of Pij, see [76]), and the minimisation equations simply read: 

mh = bi 2 tan /3 — - cos 2/3, m 2 2 = &12 cot /3 + cos 2/3, 

?n 2 3 = &32 tan/3, m 2 4 = 614 cot/3, (32) 

m i 5 = ^52 tan /3 , m 2 6 = 616 cot /3. 

After minimising the potential 7 , one can derive the charged, neutral scalar and pseudoscalar mass 
matrices, and obtain the mass eigenstates and the diagonalisation matrices. For the neutral states 

' We have verified that for each of the configurations analysed in this work, we are indeed in the presence of a local 
minimum with respect to the neutral Higgs scalars. Not only have we imposed the conditions for an extremum, but we 
also verified that it was a minimum by checking that all the minors of the Hessian matrices were positive definite. In 
terms of the Higgs spectrum, this is reflected in the absence of charged and neutral tachyonic states. Nevertheless, we 
do not discard the possibility of a global minimum associated with non-vanishing VEVs for the charged components of 
the six Higgs doublets. 
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(scalars and pseudoscalars), the relation of mass and interaction eigenstates is given by 

S r M 2 r S ] r = A 2 r = diag (mf) , i=l,...,6, 

S!M]S\ = A? = diag(mf 2 ), i = l,..y,6, (33) 

with A 2 r j the diagonal scalar and pseudoscalar squared mass eigenvalues (notice that the i = 1 term 
for the pseudoscalars corresponds to the unphysical massless would-be Goldstone boson). We recall 
here that we are working in the Higgs-basis, and that the matrices that diagonalise the mass matrices 
in the original basis can be related to the latter as 

Sa, v = S RtI P, (34) 

where P is the matrix appearing in Eq. (29). A detailed study of such an extended Higgs sector, 
including choice of basis, minimisation of the potential and derivation of the tree-level mass matrices 
can be found in [76]. It is also important to notice that throughout the analysis, and since our aim is 
to investigate to which extent FCNCs push the lower bounds on the Higgs masses, we do not consider 
radiative corrections to the Higgs masses, using the bare masses instead. 

4 Yukawa interactions and tree-level FCNCs 

In the quark and Higgs mass eigenstate basis, the Yukawa interaction Lagrangian reads: 

^Yukawa = ~ £ [ (V d )^ h* d a R d\ + i h? d% d\ + H.C.' 

v £=1,3,5 

~ [( V “)ofc h j u R u b L + i (W u )^ b h p 3 u a R u b L + H.c. . (35) 

In the above, a,b denote quark flavours, while i,j = 1,...,6 are Higgs indices, with s (p) denoting 
scalar (pseudoscalar) mass eigenstates. The latter are related to the original states as h s = S a a, 
h p = S v ip, as from Eqs. (33,34). The scalar (pseudoscalar) coupling matrices V (W) are defined as 

(V q )l =(Sih (VZY?V*U, 

=(4)« (36) 

with i = 1,3,5 (2,4,6) for q = d{u) and j = 1,...,6. Y q denote the Yukawa couplings, whose 
down-type elements were displayed in Eq. (15), and Vl,r are the unitary matrices that diagonalise 
the quark mass matrices as 

V R M q Vl ] = diag (m?), q = u,d, (37) 

so that the Cabibbo-Kobayashi-Maskawa matrix is defined as 

Vckm = V£Vf. (38) 

We emphasise that the matrices Vl, R which diagonalise the ciuark mass matrices do not, in general, 
diagonalise the corresponding Yukawa couplings. Hence, both scalar and pseudoscalar Higgs-quark- 
quark interactions may exhibit a strong non-diagonality in flavour space, which in turn translates in 
the appearance of FCNCs and CP violation at the tree-level. 
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Even though a detailed discussion of FCNCs in multi-Higgs doublet models was presented in [76], 
we summarise here some relevant points, focusing on the neutral kaon sector and investigating the 
tree-level contributions to A nriK- The latter is simply defined as the mass difference between the long- 
and short-lived kaon masses, 

A n%K = ttik l — niK s — 2 |A4^| . (39) 

The contribution to A4f 2 associated with the exchange of scalar Higgses (with masses to*) is given by 



while the exchange of a pseudoscalar state (with mass rri 1 - ) reads 



Once all the contributions to have been taken into account, the prediction of this orbifold model 
regarding A v%k should be compared with the experimental value, ( Ato^ ) eX p — 3.49 x 10“ 12 MeV [1]. 
For the other neutral meson systems, Bd , B s and D°, the computation is analogous. In each case 
the viability of the model imposes that the obtained results should be compatible with the current 
bounds: ( ArriB d )exp — 3.304 x 10“ 13 GeV, (Amsjexp > 94.8 x 10“ 13 GeV and Am. D o < 46.07 x 10“ 12 
MeV [1]. 

Before proceeding to the numerical analysis, let us briefly comment on the several contributions. 
First of all, it is widely recognised that, in models with tree-level FCNCs, the most stringent bounds 
are usually associated with A ttik- For the Z 3 orbifold scenario, with hierarchical Yukawa couplings, 
one expects the bound from A ms d to be less severe than that of A mx- The same should occur for the 
B s mass difference, since in the SM this mixing is already maximal (the only exception occurring if 
new contributions matched exactly those of the SM, but had opposite sign, in which case a cancellation 
could take place). The D° mass difference can be quite challenging to accommodate. As pointed out 
in [60,63] and [84], models allowing for FCNC at the tree-level may present the possibility of very 
large contributions to A mo, and the latter could even exceed by a factor 20 those to A wj- [63]. 
Nevertheless, one should bear in mind the fact that mixing in the D° sector is very sensitive to the 
hadronic model used to estimate the transition amplitudes, and there is still a very large uncertainty 
in deriving its decay constants, etc. Therefore, the constraints on a given model arising from A tod 
should not be over-emphasised, and we will adopt this conservative view throughout our discussion. 

Another interesting issue 8 is that of rare decays. It has been argued that, again when no theory for 
the full Yukawas is available, some rare decays become very sensitive to flavour changing contributions 

8 Extended Higgs sectors with flavour violation have other interesting consequences, such as flavour violating Higgs 
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induced by Higgs exchange at the tree-level [69]. In the present model, the Yukawas are well-defined, 
not only for the quark, but also for the lepton sector. In a forthcoming work [77], we will analyse in 
detail the lepton sector of this class of orbifold constructions, taking also into account the potentially 
most constraining decay modes, as p. —> eq, Bd —> K^lt and B s —> / it. 

Additionally, and given the existence of flavour violating neutral Higgs couplings, and the possibil¬ 
ity of having complex Yukawa couplings, it is natural to have tree-level contributions to CP violation. 
In the kaon sector, indirect CP violation is parameterised by ex, and defined as 


gjTr /4 J m | M12 A 9 

V2 |A U | 2 A m K 


(42) 


where A„ is defined from CKM elements as A„ = V* a V u d- From experiment one has ex = (2.284 ± 
0.014) x 10 -3 [1], In this case, and since we are in the presence of tree-level, rather than 1-loop 
interactions, the new contributions to ex are expected to be quite large, even if the amount of CP 
violation associated with the CKM matrix (for instance parameterised by Jcp [88 ]) is far smaller than 
the value derived from usual SM fits - Jqp ~ 0(1O -5 ) [1]. 


5 Numerical results 

In the present scenario, most of the observables addressed in the previous Section receive their domi¬ 
nant contributions from tree-level processes. This situation strongly diverges from the usual scenarios 
of both SM and MSSM, where FCNCs only occur at the 1-loop level. Given the increasing experi¬ 
mental accuracy, it is important to investigate to which extent the present scenario is compatible with 
current experimental data. 

We divide the numerical approach in two steps. Firstly, we focus on the string sector of the model, 
and for each point in the space generated by the free parameters of the orbifold (£ 5 , a/), we derive 
the up- and down-quark mass matrices 9 and compute the CKM matrix. This procedure allows us 
to investigate the several regimes of parameters that translate into viable quark spectra, and discuss 
the implications of the relations between the several parameters. At this early stage, we consider 
only real values for the orbifold parameters. Further imposing the conditions associated with EWSB 
given in Eq. (18), and fixing a value 10 for tan/3, one can then determine the values of g N and £1 (cf. 
Eqs. (20,21)). Another possible approach would be to scan over the space generated by the moduli (7/) 
and the VEVs of the 5/7(3) x 5/7(2) singlet fields (d), but this would translate in less straightforward 
relations between the orbifold parameters and the experimental data. A secondary step requires 
specifying the several Higgs parameters, which must obey the minimum criteria of Eq. (32). Finally, 
the last step comprehends the analysis of how each of the Yukawa patterns constrains the Higgs 
parameters in order to have compatibility with the FCNC data. In particular, we want to investigate 

and top-decays. For a discussion of the latter, and the associated experimental signatures at the next generation of 
colliders, see, for example [85-87], and references therein. 

9 It is worth emphasising here that the quark masses appearing in Eq. (16) (and in all subsequent relations) are used 
on the sole purpose of obtaining an approximate determination of the VEVs. Afterwards, the actual values of m qi are 
exactly computed. 

10 We recall here that tan (3 is a necessary parameter to specify the Higgs sector, which in turn is mandatory to 
investigate the issue of FCNCs in the present model. 
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Set 

m u 

m d 

m c 

m s 

m t 

mb 

A 

0.004 

0.008 

1.35 

0.13 

180 

4.4 

B 

0.0035 

0.008 

1.25 

0.1 

178 

4.5 

C 

0.0035 

0.004 

1.15 

0.08 

176 

4.1 

D 

0.004 

0.006 

1.2 

0.105 

178 

4.25 


Table 1: Sets of input quark masses (in GeV) used in the numerical analysis. 

how heavy the scalar and pseudoscalar eigenstates are required to be in order to accommodate the 
observed values of A uik, A m,B d , etc. 

5.1 Quark Yukawa couplings and the CKM matrix 

As discussed in [25], there are three regimes for the values of cd and (3 *, depending on the specific 
orbifold configuration : 

(a) cd ~ £5 and ~ 1 ; 

(b) a? w £5 and ss l/£ 5 ; 

(c) a/ ~ 1 and fif ~ 1 . 

In any case, it is clear from Eq. ( 8 ) that oJ and fif must obey, by construction, the following bounds: 

£5 < a f < 1, 1 < P f < — ■ (43) 

£5 

In what follows we investigate whether each point in the orbifold parameter space can be associated 
with a consistent quark spectrum and mixings. For given values of the input quark masses, one fixes 
the ratio of the several Higgs VEVs, which in turn allows to reconstruct the full quark mass matrices, 
and obtain the mass eigenstates and CKM matrix. In particular, throughout this analysis we shall 
focus on four sets of input quark masses, whose values are listed in Table 1. 


Throughout, we require the CKM matrix elements to lie within the following ranges [1]: 

/ 0.9739- 0.9751 0.221 - 0.227 0.0029 - 0.0045 \ 

Vckm = 0.221 - 0.227 0.9730- 0.9744 0.039-0.044 . (44) 

v 0.0048- 0.014 0.037- 0.043 0.9990-0.9992 j 

Regarding the effect of the orbifold parameters on quark mixing, let us notice that both £5 and a u are 
crucial in obtaining the Cabibbo angle. As expected, the down-sector parameters are those directly 
responsible for the mixings between generations, and their role is particularly relevant in determining 
Vtd ( a dC - and to a lesser extent, also a M<= ), Vt s , V u b and V c b- Once these elements are fixed in 
accordance with experiment, the others are usually also in agreement. Finally, let us recall that 
from choosing a concrete value for tan /?, and complying with the bound on Mz from EW symmetry 
breaking, Eqs. (18,20), one is implicitly fixing for each set of {£ 5 , a/,/?Q, the values of £1 and g N. 
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Figure 2: Correlation between the orbifold parameters (a“°, £ 5 ) and (a dC , £ 5 ), for distinct sets of input 
quark masses, A-D (red full lines, red dashed lines, full black and dashed black lines, respectively). 


In Fig. 2, we present the correlation between the orbifold parameters, for the four sets of input 
quark masses given in Table 1. We only present points that are associated with viable quark masses 
and that are in agreement with current bounds on |Vckm| (from Eq. (44)). Rather than scatter plots, 
in Fig. 2 we present sets of points. This is due to having very narrow intervals of fluctuation for all 
the parameters. As an example, let us mention that for constant values of £ 5 , a d and (3^, a u is 
fixed within a 2% interval. ^From Fig. 2 it is clear that for a given set of Higgs VEVs (determined 
by the associated set of input quark masses - Table 1) the allowed orbifold parameter space is very 
constrained. This is a direct consequence of the increasing accuracy in the experimental determination 
of the Vckm- In each set (A D), moving outside the displayed ranges would translate in violating 
the experimental bounds on (at least) one of the CKM matrix elements. Larger values of £5 would 
also (typically) be associated with an up-quark mass below the current accepted range. Regarding the 
values of (3*, these are constrained to be (3^ ss 1 throughout the allowed parameter space (cf. Eq. (12)). 
From the direct inspection of Fig. 2, together with the fact that (3? ~ 1, it appears that at least two 
regimes for af are present. For the up sector, one would suggest that the orbifold configuration is 
such that is roughly 0[e 5 ), so that one is faced with regime (a). The same would happen for sets 
A and B in the down sector, although sets C and D appear to favour a regime with ~ 1, thus 
suggesting regime (c). 

Let us now aim at understanding the behaviour of both a u ,d as a function of £ 5 . In Ref. [25], 
several analytical relations for the CKM matrix elements as a function of {£ 5 , 0 /,/?!} were derived. 
Although those relations were computed for the case of hermitic mass matrices, where V[ = V^, 
and are thus not truly valid for the present case, they are quite useful in understanding Fig. 2. For 
instance, the Cabibbo angle is given, to a very good approximation, by 


V u . 


2 ( mt 1 
e 5 I- 77c 

\m c a u 


m b 1 \ 

H c I ’ 
m s or J 


(45) 


and the above expression is clearly dominated by the first term on the right-hand side (r.h.s). From 
Eq. (45), it becomes transparent that the dependence of a u on £5 should indeed be parabolic, as 
clearly displayed in Fig. 2. Regarding a d , its evolution is strongly influenced by the allowed regions 


15 




















0.01 0.015 0.02 0.025 



0.25 

0.2 

0.15 

0.1 

0.05 

0 


Figure 3: Values of a d and a u (full and dashed lines, respectively) as a function of £ 5 . Vertical 
dotted lines isolate the ranges of £5 associated with each set A-D. 


of £5 (as dictated by the quark mass regimes taken as input, especially m s /m.b). For sets A and B, 
the ratio of down-type quark masses is such that the leading term in the analytical expression of V u b, 

„d c 


V ub 


^a dC - —a" 
m b m t 


- £5 


m d ni t or 
mb m c a uC (3 dc 


m u 1 
m t /3“‘ 


o(4), 


(46) 


provides a reasonable understanding of a d (£ 5 ). For sets C and D, the situation is more involved, and 
the second term on the r.h.s. of Eq. (46) plays an important role. In fact, a d receives contributions 
which display a near-resonant behaviour for the input quark mass ratios in the correspondent £5 range. 
This is the origin of the unexpected positioning of set C in Fig. 2. 

We now address the behaviour of a u and a d . By construction, and as clearly seen from Eq. (11), 
once £5 and oJ (or equivalently £5 and fjf) are set, one is implicitly fixing aV Moreover, satisfying 
the EWSB conditions for the VEVs, together with imposing a given value of tan/3 also translates in 
determining g N and e\. In Fig. 3, we plot the values of a u and a d as a function of £ 5 , as determined 
from Eq. (11). As seen from Fig. 3, the behaviour of set C regarding a d is slightly misaligned with 
what one would expect from the analysis of sets A, B and D. Notice however that this is due to the 
dependence of a d on a d (cf. Eq. (11)). Recall from Fig. 2 that for set C, the allowed values of a d 
were somewhat higher than for the other cases, and this is the source of the suppression displayed in 
Fig. 3, set C. 


In Fig. 4, we present the values of £1 as a function of £5 for the four sets of quark masses, and 
distinct regimes of tan/3. ^From Fig. 4 we can also verify that the values of £1 are, in general, larger 
than those of £ 5 . The “misaligned” behaviour of set C is again a consequence of the effects already 
discussed. Regarding the actual value of gN it suffices to mention that for the orbifold parameter 
space here investigated, and for the values of tan/3 here selected, one has 1.03 < g N < 1.16. 

Since we now have the relevant information, we can further compute the value of the lattice’s 
diagonal moduli, T) and T 5 , as defined in Eq. (4). The value of T 5 is unambiguously determined. 
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Figure 4: e\ as a function of £5 for tan /? =3, 5, 10 and 20 (red full lines, red dashed lines, full black 
and dashed black lines, respectively). Vertical dotted lines isolate the ranges of £5 associated with 
each set A-D. 
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Figure 5: Diagonal lattice moduli, Xi (left) and X 5 (right), as a function of £ 5 . For the case of Ti, we 
plot it for several values of tan/? =3, 5, 10 and 20 (following the line and colour scheme of Fig. 4). 
For X 5 the dotted line denotes the prediction of the orbifold. 


Nevertheless, and since the determination of £1 is a direct consequence of complying with the EWSB 
conditions for a fixed value of tan/?, its determination varies accordingly. In Fig. 5 we display the 
diagonal moduli as a function of £ 5 , for tan/? =3, 5, 10 and 20. ^From Fig. 5, it is manifest that for 
the parameter space investigated, the values of the diagonal moduli, Xi and X 5 , are never degenerate. 
This confirms our original assumption (see footnote 2) that distinct moduli are indeed required to 
accommodate the experimental data. Although we have allowed for non-degenerate X^, this remains 
quite a restrictive choice. We recall that we still have six additional degrees of freedom, namely the 
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Figure 6 : Down- and up-type Higgs VEVs, w 1 , 3,5 and W 2 , 4 , 6 > as a function of £ 5 , for tan/3 =5. For 
convenience, we plot we/10- As before, vertical dotted lines isolate the ranges of £5 associated with 
each set A-D. 


angles between the complex planes, which we have taken as zero in the present analysis. 

To conclude the study of the orbifold parameter space, let us just plot the values of the Higgs 
VEVs, again as a function of £ 5 . As an illustrative example, depicted in Fig. 6 , we take tan/3 = 5. It 
is interesting to comment that in the up-sector, the VEVs exhibit a clearly hierarchical pattern, W 2 < 
W4 < wg while in the down-type VEVs we encounter a not so definite pattern, with W3 < W5 < w±. 
This is a direct consequence of the relations given in Eq. (16). 

Finally, let us summarise our analysis of the orbifold parameter space by commenting on the 
relative number of input parameters and number of observables fitted. Working with the six Higgs 
VEVs (wi), and the orbifold parameters £ 1 , £ 5 , and one can obtain the correct EWSB ( Mz ), 
as well as the correct quark masses and mixings (six masses and three mixing angles). 

5.2 Tree-level FCNCs in neutral mesons 

The present orbifold model does not include a specific prediction regarding the Higgs sector. For 
instance, we have no hint regarding the value of the several bilinear terms, nor towards their origin. 
Concerning the soft breaking terms, the situation is identical. Since the FI D-term, which could have 
broken SUSY at the string scale, was cancelled, one must call upon some other mechanism to ensure 
that supersymmetry is indeed broken in the low-energy theory. In the absence of further information, 
we merely assume that the structure of the soft breaking terms is as in Eq. (26), taking the Higgs soft 
breaking masses and the B^i -terms as free parameters (provided that the EW symmetry breaking and 
minimisation conditions are verified). 

Before proceeding, it is important to stress that the Higgs spectrum (i.e. the scalar and pseu¬ 
doscalar physical masses) cannot be a direct input when investigating the occurrence of FCNCs. In 
some previous studies of FCNCs in multi-Higgs doublets models (see for example [63]) the bounds 
were derived for the diagonal entries in the scalar and pseudoscalar mass matrices. However, this 
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approach neglects mixings between the several fields Ri , It , and excludes the exchange of some scalar 
and pseudoscalar states. 

Although it is possible to begin the analysis from the original basis (where all neutral Higgses 
develop VEVs), we prefer to define the Higgs parameters on the Higgs-basis, relying on the min¬ 
ima conditions (and associated inequalities) to ensure that we are on the presence of true minima. 
Therefore, the parameters that must be specified are: 

tan/3, rn-j, fey, (47) 

entering in Eq. (32). In the absence of orbifold predictions for the Higgs sector parameters, and 
motivated by an argument of simplicity, we begin our analysis by considering textures for the above 
parameters as simple as possible. 

To avoid rewriting the Higgs soft-breaking masses, we adopt the following parameterisation: 


( ® 

<g> 

« \ 

( ® 

0 

0 \ 


X3 

y xlTeV, = 

® 

X4 

y x ITeV, ^/b~ = b x ITeV. (48) 

^ <s> 

y 

) 

^ <8> 

y 

x& ) 


In the above, m d( ' u > should be understood as the i.j = 1,3,5 (k, l = 2,4,6) submatrices of the 6 x 6 
matrix that encodes the rotated soft-breaking Higgs masses in the Higgs basis (see [76]). The symbol 
0 denotes an entry which is fixed by the minima equations of Eq. (32). This parametrisation allows to 
easily define the Higgs sector via six dimensionless parameters. We begin by taking a near-universality 
limit for the Higgs-sector textures introduced in Eq. (48). Regarding the value of tan/3, and if not 
otherwise stated, we shall take tan /3 = 5 in the subsequent analysis. We first consider the following 
four cases, presenting the associated scalar and pseudoscalar Higgs spectra: 

(a) x$ — X 4 — 0.5 , X 5 = xq = 0.75, y = 0.1, b = 0.1 
m s = {82.5,190.6,493.9,515.9,744.4,760.2} GeV; 
m p = {186.8,493.9,515.9,744.4,760.2} GeV. 

(b) £3 = £4 = 0.4 , £5 = Xq = 0.8, y = 0.15 , b = 0.2 
m s = {84.6,213.9,387.4,560.8,785.9,879.1} GeV; 
m p = {215.2,387.3,560.5,785.9,878.9} GeV. 

(c) X 3 = X 4 = 0.75 , X 5 = xq = 1, y = 0.25, b = 0.2 

m s = {83.6,292.9,733.6,785.9,987.6,1057.0} GeV; 
m p = {291.1,733.6, 785.9,987.6,1057.0} GeV. 

(d) X 3 = X 4 = 0.3 , £5 = xq = 0.8, y = 0.1, b = 0.1 

m s = {79.4,121.5,296.9,354.3,794.6,808.8} GeV; 
m p = {114.8,296.9,353.7,794.6,808.8} GeV. 

In Fig. 7 we plot the ratio Auik /(Amj{)exp versus £- 5 , for cases (a)-(d). We considered tan/3 = 5 
and, since the other sets of input quark masses were in general associated to smaller FCNC effects, 
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Figure 7: Am*-/(Amif)exp as a function of £5 for set B and tan (3 = 5. The Higgs parameters 
correspond to textures (a)-(d). 

we take the quarks masses as in “set B” (Table 1). Once again, all the points displayed comply 
with the bounds from the CKM matrix. /.From Fig. 7 it is clear that it is quite easy for the orbifold 
model to accommodate the current experimental values for Am*-. Even though the model presents the 
possibility of important tree-level contributions to the kaon mass difference, all the textures considered 
give rise to contributions very close to the experimental value. Although (a) and (b) are not in 
agreement with the measured value of A rriK, their contribution is within order of magnitude of 
(Amic)exp- As seen from Fig. 7, with a considerably light Higgs spectrum (i.e. m h 0 < 1 TeV), one is 
safely below the experimental bound, as exhibited by cases (c) and (d). This is not entirely unexpected 
given the strongly hierarchical structure of the Yukawa couplings (notice from Eq. (15) that Y^i is 
suppressed by £§). 

One important aspect clearly manifest in Fig. 7, and which has been overlooked in some previous 
analyses, is that the Higgs mixings can be as relevant as the Higgs eigenvalues in determining the 
contributions to A uik- This is patent in the comparison of curves (c) and (d). From a naive inspection 
of the Higgs spectra associated to each case, one would expect that (c) would induce a much stronger 
suppression to the model’s contribution to ISjuk- Nevertheless, case (d), with a spectrum quite 
similar to case (b), and indeed much lighter than that of (d), but with much smaller mixings, is the 
one associated with the strongest suppression of A mx- It is also important to comment on the effect 
of changing tan/3 regarding the contributions to the kaon mass difference. For the specific case of 
texture (c), let us investigate the effect of varying tan/3. We take the quark masses as in set B, and 
while keeping the Higgs parameters fixed, tan/3 is taken in the range 3 < tan/3 < 9.5. As it becomes 
clear from Fig. 8 , larger values of tan/3 produce increasingly larger contributions to A mx- This is a 
direct consequence of the fact that, due to larger off-diagonal terms in the Higgs mass matrices, the 
mixing is larger, and the corresponding eigenstates become lighter. Even though the masses of the 
heaviest states remain stable, the intermediate states become lighter, and the FCNC contributions 
are less suppressed. Close to tan/3=10, it is no longer possible to find physical minima of the Higgs 
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Figure 8: ArriK / (ArriK)exp as a function of tan/3 for set B. The Higgs parameters correspond to 
texture (c). 


potential, and tachyonic states emerge. This effect has been already pointed out in the general analysis 
of [76]. 

5.2.1 Bd and B s meson mass difference 

For the several parameterisations of the Higgs sector used for the analysis of A niK, we display in 
Fig. 9 the contributions of Higgs textures (a)-(d) to the B d mass difference. As one would expect, 
given the structure of the mass matrices (and thus the Yukawa couplings), the present model generates 
very small contributions to ArriB d - All the textures analysed, even those associated with excessive 
contributions to A?r Ik are in good agreement with the experimental data on the B d mass difference. 
Notice that the behaviour of the textures is now quite distinct: as an example, texture (b), which 
generated the second largest contribution to the A rriK is now the one associated with the strongest 
suppression. This stems from the fact that the leading contributions are now given by distinct Higgses, 
whose couplings to the quarks need not be identical. 

Likewise, in Fig. 10 we plot the contributions to the B s mass difference. In this case experiment 
only provides a lower bound, so that any ratio larger than 1 is in agreement with current data. As 
we would expect from the discussion of Section 4, this model’s contributions to A ms s are well above 
the current limit. 


5.2.2 D° — D° mass difference 

The cases (a)-(d) considered in the previous subsections generate contributions to A mo that exceed 
the experimental bounds by at least a factor 10. As discussed in Section 4, this is not surprising, nor 
excessively worrying. Nevertheless, and for the sake of completing the analysis, let us consider four 
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Figure 9: Amj 3 d /(Am Bd ) exp as a function of £5 for set B and tan/3 = 5. 
correspond to textures (a)-(d). 
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Figure 10: Am_B s /(Am_B s ) e xp as a function of £5 for set B and tan/3 = 5. The Higgs parameters 
correspond to textures (a)-(d). 


additional Higgs patterns, in order to derive a bound on the mass of the heaviest Higgs boson that 
would render the model compatible with the data from the D° sector. 

The new Higgs textures are defined as follows: 

(e) £3 = X 4 = 0.75 , £5 = 2.5, xe = 7.5 , y = 0.5 , b = 0.5 
m s = {84.2,672.7,704.9,1414,2573,7501} GeV; 
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Figure 11: Amo/{Am,£>) exp as a function of £5 for set B and tan/3 = 5. The Higgs parameters 
correspond to textures (e)-(h). 

m p = {673.1, 704.9,1414,2573,7501} GeV. 

(f) X 3 = X 4 = 0.5 , X 5 = 5, xq = 7.5, y = 0.5 , b = 0.1 
m s = {82.7,201.4,492.4,516.4,5000, 7500} GeV; 
m p = {197.9,492.4,516.4,5000, 7500} GeV. 

(g) X 3 = x± = 1, X 5 = X 6 = 7.5, y = 0.5 , b = 0.5 
m s = {82.3,378.5,958.9,1578, 7483,7518} GeV; 
m p = {379.2,958.9,1578, 7484,7518} GeV. 

(h) X 3 = x 4 = 3 , x 5 = x 6 = 7 , y = 1, b = 0.5 

m s = {84.1,1032,2964,3059,6984,7022} GeV; 
m p = {1031,2964,3059,6984,7022} GeV. 

For the cases (e) to (h) we present in Fig. 11 the contributions to Amo- Compatibility with experiment 
can be obtained for any of the sets (e), (g) or (h), thus suggesting that one of the Higgses (a state 
mostly dominated by an up-type Higgs field) must be at least of around 7.5 TeV. Notice that no 
major hierarchy is required from the Higgs spectrum - case (e) is an excellent example of the latter 
statement, in the sense that one obtains states softly ranging from 600 to 7500 GeV. One may wonder 
if such a choice of Higgs soft-breaking terms may lead to a fine-tuning problem. In [76], it was pointed 
out that for non-degenerate VEVs, soft-breaking terms above the few TeV range typically induced a 
fine tuning stronger than 1%. Nevertheless, we again stress that these values for the Higgs masses (as 
derived from the D° mass difference analysis) should not be interpreted from a very strict point of 
view. 
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Figure 12: £/f/(£/f)exp as a function of |£s| for set B and tan/3 = 5. The Higgs parameters correspond 
to textures (c)-(e). 


5.3 Tree-level CP violation: ek 

Finally, we turn our attention to the issue of CP violation. So far, we have seen that accommodating 
the several A F = 2 observables is not an excessively hard task (especially if we choose to set aside the 
D° sector). Nevertheless, a successful model of particle physics must also comply with the observed 
CP violation in the kaon sector. As we mentioned in Section 4, we will only consider the specific 
contribution of the present model to indirect CP violation in the kaon sector ( ek )• Therefore, we now 
assume £5 (and thus a?) to be a complex quantity, and parameterise it as £5 = |es| e 1 ^. 

In Fig. 12 we present the tree-level contributions to ek (normalised by its experimental value) 
as a function of |£s|. We take the input quark masses as in set B, fix tan/3 = 5, and analyse the 
Higgs textures associated with cases (c), (d) and (e). For each texture, the phase is assumed to be 
4> c = 4.0 x 10~ 4 , 4>d = 2.5 x 10 -4 and <j> e = 2.5 x 10 -3 . These phases are taken as illustrative examples; 
we choose values that for a specific set of input quark masses (set B, in this case) and a given Higgs 
texture (c)-(e) simultaneously succeed in generating an amount of ek close to the value experimentally 
measured (range delimited by dotted grey lines in Fig. 12), and still have a compatible CKM. One 
should bare in mind that once arges (and thus arg 0 ^) is no longer a small number, it will significantly 
affect the computation of the Yukawa couplings, and thus the CKM matrix. We have considered one 
texture that accommodates all FCNC observables, namely texture (e), which already has a somewhat 
heavy Higgs spectrum. In this case, the phase required to saturate (eif)exp is 0(1O~ 3 ). Since we 
do not wish to view the D° sector as a very stringent constraint, we also consider two other Higgs 
patterns, (c) and (d), which only succeed in complying with both kaon and B-meson data. In these 
cases, the phases necessary to obtain (£_K-) exp are now (D( 10 -4 ), as one would expect, since a lighter 
set of Higgs bosons typically enhances the contributions. 

In the range of parameters analysed in this plot, the amount of CP violation stemming from the 


24 








CKM matrix is Jcp ~ 0{ 10 ~ 8 — 10~ 6 ), i.e. at least one order of magnitude below the SM value 
that is associated with the observed Sk [1]- The possibility of obtaining an orbifold configuration 
that saturates the observed value of £/<- and at the same time allows to reproduce the Jcp required 
by the unitarity triangle fits should not be discarded. It is clear that the phase of £5 must be quite 
large, and such values would push us to distinct areas of the orbifold parameter space. It is worth 
emphasising that there are still other sources of CP violation in addition to the one we have studied 
in this Section. As mentioned in footnote 3, a more complicated choice of the VEVs c,; could lead to 
physical phases in the quark masses matrices, which would in turn contribute to the physical CKM 
phase. 


6 Conclusions 

In this work we have investigated whether it is possible to find abelian Z 3 orbifold configurations that 
are associated with an experimentally viable low-energy scenario. 

This class of models provides a beautiful geometric mechanism for the generation of the fermion 
mass hierarchy. The Yukawa couplings are explicitly calculable, and thus a solution to the elusive 
flavour problem of the SM and MSSM can in principle be obtained. 

We have concentrated here in Z-$ orbifold compactifications with two Wilson lines, which naturally 
include three families for fermions and Higgses. The fact that additional Yukawas are present opens 
the possibility of obtaining realistic fermion masses and mixings, entirely at the renormalisable level 
(with a key role being played by the FI breaking). We have surveyed the parameter space generated 
by the free orbifold parameters, and we have verified that one can find configurations that obey the 
EW symmetry breaking conditions, and can account for the correct quark masses and mixings. 

On the other hand, having six Higgs doublets (and thus six quark Yukawa couplings) poses the 
potential problem of having tree-level FCNCs. By assuming simple textures for the Higgs free param¬ 
eters, we have verified that the experimental data on the neutral kaon mass difference, as well as on 
A vriB d and A mp s can be easily accommodated for a quite light Higgs spectra, namely m h 0 < 1 TeV. 
The data from the D° sector proves to be a more difficult challenge, requiring a Higgs spectrum of 
at least 7 TeV, but we again stress that, in view of the theoretical and experimental uncertainties 
associated with the D° sector, this constraint should not be over-emphasised. 

CP violation can be also embedded into the low-energy theory. Although CP is a gauge symmetry 
of the full theory, it can be spontaneously broken at the string scale, if the VEVs of the moduli have a 
non-vanishing phase. We have parameterised these effects by assuming the presence of a phase in £ 5 , 
and we have verified that one can also obtain a value for Sk in agreement with current experimental 
data. 

The presence of a fairly light Higgs spectrum, composed by a total of 21 physical states, may 
provide abundant experimental signatures at future colliders, like the Tevatron or the LHC. In fact, 
flavour violating decays of the form hi —> qq, or hi —* l + l~ may provide the first clear evidence 
of this class of models. Z 3 orbifold compactifications with two Wilson lines are equally predictive 
regarding the lepton sector. This analysis, especially that of the neutrino sector, will be addressed in 
a forthcoming work [77]. 
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